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The  dissociative  recombination  of  Na+(D20)  ion  has  been  studied  at  the  heavy-ion  storage  ring 
CRYRING  (Manne  Siegbahn  Laboratory,  Stockholm  University).  The  cross  section  has  been 
measured  as  a  function  of  center-of-mass  energy  ranging  from  1  meV  to  0. 1  eV  and  found  to  have 
an  E~131  dependence.  The  rate  coefficient  has  been  deduced  to  be  (2.3±0.32) 
X  10~7(re/300r°-95±°-01  cm3  s"1  for  Te—  50- 1000 K.  The  branching  ratios  have  been  measured  at 
0  eV.  Of  the  four  energetically  accessible  dissociation  channels,  three  channels  are  found  to  occur 
although  the  channel  that  breaks  the  weak  Na+-D20  bond  is  by  far  dominant.  ©  2004  American 
Institute  of  Physics.  [DOI:  10.1063/1.1812532] 


I.  INTRODUCTION 


Na+  (H20)  may  also  result  from  raptcnigana  swncnmg  from 


Easily  ionizable  alkali  metals  emitted  from  surfaces  can 
contribute  significantly  to  the  plasma  that  forms  around  hy¬ 
personic  (>Mach  5)  vehicles  during  atmospheric  reentry. 
Due  to  the  high  pressure  of  atmospheric  gases  around  the 
vehicle  and  the  low  radiative  recombination  rate  of  alkali 
atomic  ions,  plasma  neutralization  is  dominated  by  three- 
body  processes,  e.g., 


Na+(N2)  or  Na+(C02).  This  is  certainly  the  case  in  the 
atmosphere  for  H30+(H20)  formation  from  H30+(N2).16 

The  second  step  in  the  chaperone  mechanism  involves 
the  ionic  cluster  recombining  with  a  free  electron, 


Na+(H20)  +  e  — > neutral  products.  (3) 


Na+  +  e~  +  M— >Na+M,  (1) 

where  M  is  an  atmospheric  gas  or  an  electron.  Alternatively, 
a  third  body  can  increase  the  recombination  rate  through  a 
chaperone  mechanism,  which  takes  place  in  two  steps.  The 
alkali  ion  first  clusters  to  an  atmospheric  gas  molecule  such 
as  water, 

Na+  +  H2  O  +  M<-»  Na+  ( H2  O)  +  M.  (2) 

Although  it  might  be  instructive  to  consider  clustering  to  the 
most  abundant  atmospheric  species,  N2  and  02,  Na+(H20) 
is  expected  to  be  more  abundant  than  the  Na+  •  N2  and 
Na+  •  02  complexes  due  to  much  stronger  binding  energy  in 
Na+(H20).1,2  The  interaction  between  Na+  and  water  has 
been  the  subject  of  numerous  experimental3,4  and 
theoretical5-14  studies  due  to  its  important  role  in  many 
chemical  and  biochemical  phenomena.  The  formation  of 
Na+(H20)  in  the  presence  of  a  third  body  [reaction  (2)] 
proceeds  in  two-step  collision  processes,  where  the  first  step 
is  water  attachment  to  the  ion,  after  which  the  cluster  is  left 
unstable,  and  the  second  step  is  removing  excess  energy  by 
means  of  collision  with  the  third  body,15  i.e.,  the  well  known 
Lindeman  mechanism  for  association.  The  formation  of 


This  is  called  dissociative  recombination  (DR)  and  is  usually 
much  faster  than  the  three-body  process  shown  in  reaction 
(1).  The  equilibrium  constant  for  reaction  (2)  is  known3  but 
the  rate  for  DR  of  Na+(H20)  [reaction  (3)]  is  not.  Here,  we 
report  a  study  of  reaction  (3)  made  at  the  heavy-ion  storage 
ring  CRYRING  to  better  assess  the  importance  of  the  chap¬ 
erone  mechanism.  In  order  to  better  obtain  product  distribu¬ 
tions,  Na+(D20)  was  used  to  compensate  for  the  low  reso¬ 
lution  of  the  detector.  The  isotopic  shifts  in  the  geometry, 
rate  constants  and  branching  ratios  are  expected  to  be  small 
based  on  previous  experience  with  systems  of  similar  size. 

The  dissociative  recombination  of  Na+(HzO)  is  also  rel¬ 
evant  to  the  chemistry  of  the  D-region  of  the  ionosphere, 
where  a  10  km  half-width  Na+  belt  is  present  centered  near 
93  km  altitude,17  and  to  the  formation  of  sporadic  sodium 
layers,  which  can  occur  at  the  altitudes  between  90  and  105 
km.18,19  Moreover,  positive  ions  of  mass  41  amu,  presumably 
Na+(H20),  have  been  observed  with  a  rocket-borne  mass 
spectrometer  at  an  altitude  of  82  km.20  Similar  consider¬ 
ations  as  discussed  above  apply  to  this  situation,  i.e.,  Na+ 
cannot  recombine  but  the  water  cluster  would. 


0021-9606/2004/121  (21  )/1 0483/6/S22.00 


10483 


«sr*  #"»** 


©  2004  American  Institute  of  Physics 

c 


A K  Wmi  OPktXA  fr/\  4  AfL  -Id'l  AAA  HfSmr*  emhiaH1  */\  AID  i^am 


■50484  J.  Chem.  Phys.,  Vol.  121,  No.  21,  1  December  2004 


FIG.  1 .  Layout  of  the  heavy-ion  storage  ring  CRYRING  at  Manne  Siegbahn 
Laboratory,  Stockholm. 

II.  EXPERIMENT 

Experiments  were  performed  at  the  heavy-ion  storage 
ring  CRYRING  located  at  Manne  Siegbahn  Laboratory, 
Stockholm  University,  Sweden.  A  detailed  description  of  the 
experimental  setup  can  be  found  elsewhere.21  The  layout  of 
CRYRING,  which  is  51.6  m  in  circumference,  is  shown  in 
Fig.  1.  The  ions  were  produced  in  a  hot  filament  ion  source 
(MINIS)  by  mixing  Na  from  evaporated  NaCl  with  D20: 

Na+  +  D20— >Na+(D20)+  1.04  eV.  (4) 

Deuterated  NaOD  and  Na+(D20)  are  known  to  have  the 
same  structure,22  for  which  the  geometry  is  shown  in  Fig. 
2. 14  Once  formed  in  the  source,  the  ions  were  extracted  at  40 
keV,  mass  selected  and  injected  into  the  ring  where  the  beam 
was  further  accelerated  by  a  radio  frequency  (rf)  cavity  to  the 
maximum  beam  energy,  £&  =  2.19  MeV,  which  was  limited 
by  the  magnetic  rigidity  of  the  storage  ring.23  The  data  ac¬ 
quisition  was  started  1  s  after  the  ion  injection,  which  en¬ 
sured  infrared  vibrational  relaxation  of  the  ions  to  a  thermal 
distribution  of  vibrational  states.  The  interaction  between 
ions  and  a  high-density  electron  beam  having  a  current  of  0.4 
mA  occurred  over  a  length  of  0.85  m  in  the  electron  cooler. 
The  electron  velocities  in  the  center-of-mass  frame  are  char¬ 
acterized  by  a  flattened  Maxwellian  distribution  with  parallel 
and  transverse  temperatures  of  lcTn  =  Q.l  meV  and  kTL 
=  2  meV,  respectively.24  The  neutral  particles  arising  from 
dissociative  recombination  and  collisions  with  rest  gas  mol- 


Zhaunerchyk  et  al. 

ecules  were  separated  from  the  charged  particles  in  the  di¬ 
pole  magnet  after  the  electron  cooler,  and  were  detected  by 
an  ion-implanted  surface  barrier  detector  (SBD)  with  3000 
mm2  active  area  located  at  a  distance  of  3.85  m  from  the 
midpoint  of  the  electron  cooler.  The  output  signal  from  the 
detector  was  amplified,  shaped,  and  recorded  by  a  multichan¬ 
nel  analyzer  (MCA)  giving  a  pulse-height  spectrum  showing 
the  number  of  counts  versus  the  energy  of  the  incoming  par¬ 
ticle.  Since  all  fragments  travel  with  the  ion  beam  velocity, 
the  energy  of  the  individual  neutral  particles,  En  can  be  ex¬ 
pressed  as 

E„=m„/m  Na+(D2o)x^A>  (5) 

where  m„  is  the  mass  of  the  fragment  and  »iNa+(D,o)  's  the 
mass  of  the  parent  ion.  Particles  from  the  same  recombina¬ 
tion  event  hit  the  detector  simultaneously  for  our  purposes. 

The  pressure  in  the  storage  ring  was  low 
(~lCT11Torr)  which  ensured  a  small  number  of  collisions 
with  the  background  gas  and  resulted  in  a  lifetime  of  the  ion 
beam  of  5.6  s.  This  background  was  measured  and  the  data 
corrected  accordingly.  The  cycle  time  was  chosen  to  be  7.7  s, 
which  included  ion  injection  and  acceleration  (1  s),  vibra¬ 
tional  relaxation  of  the  ions  (1  s),  data  acquisition  (3  s),  and 
resetting  of  the  ring  (2.7  s). 

III.  DATA  ANALYSIS  PROCEDURE 

A.  Dissociative  recombination  cross  section  and  rate 
coefficient  measurements 

The  collision  energy  in  the  center-of-mass  frame  can  be 
defined  as 

^Ecm  E cat  h  'JE  cooh  (6) 

where  q  is  the  elementary  charge,  Ucalh  is  the  effective  cath¬ 
ode  voltage  at  which  the  electrons  are  produced,  and  Ecool  is 
the  electron  beam  energy  when  the  average  velocities  of  two 
beams  are  matched,  which  quantitatively  equals  the  value 
(w(,/wNa+(D^0))£'4=31.3  eV,  where  me  is  the  electron  mass. 
We  have  denoted  the  subscript  “cool,”  because  at  these  con¬ 
ditions  for  lighter  ions  a  cooling  process  is  observed  result¬ 
ing  in  quenching  of  the  energy  spread  of  the  ions  by  means 
of  Coulomb  interactions.  For  detailed  description  of  the  elec¬ 
tron  cooling  see  Ref.  24.  For  an  ion  as  heavy  as  Na+(D20), 
cooling  is  not  very  effective. 

The  measurement  for  the  DR  cross  section  versus  the 
interaction  energy  was  performed  by  linearly  ramping  the 
cathode  voltage  from  43.5  to  21.2  V,  covering  the  range  of 
collision  energy  between  0  and  1  eV  twice  [Fig.  3(a)].  In 
order  to  count  the  number  of  DR  events,  the  output  of  the 
SBD  was  connected  to  a  multichannel  scaler  (MCS),  which 
recorded  the  number  of  pulses  as  a  function  of  time  [Fig. 
3(b)],  The  DR  rate  coefficient  was  detennined  from 

dN  vei ;.e2  nr] 

n _  e  e  /«\ 
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FIG.  3.  (a)  Electron  cooler  cathode 
voltage  vs  time,  (b)  The  corresponding 
DR  count  rate. 


where  dNIdt  is  the  count  rate,  ve  and  are  the  mean  elec¬ 
tron  and  ion  velocities,  Ie  and  /,  are  the  electron  and  ion 
beam  currents,  re  is  the  electron  beam  radius,  and  /  is  the 
electron  cooler  length. 

Since  a  value  of  the  measured  DR  rate  coefficient  is  the 
velocity  weighted  cross  section  averaged  over  the  electron 
velocity  distribution,25  a  Fourier  transformation  procedure  is 
applied  to  extract  the  cross  section.21  This  method  is  efficient 
for  low  collision  energy,  because  it  takes  into  account  the 
velocity  spread  of  electrons.  However,  for  the  interaction  en¬ 
ergies  in  this  study  the  velocity  spread  is  negligible,  which 
makes  it  appropriate  to  use  the  following  formula  to  calcu¬ 
late  the  cross  section: 

R 

(8) 

.  v  cm 

where  vcm  =  \j2Ecm/me  is  the  electron  velocity  in  the  center- 
of-mass  frame.  The  data  were  finally  corrected  for  the  space 
charge  of  the  electron  beam26  and  the  effect  of  the  toroidal 
magnetic  field  in  the  curved  edges  of  the  electron  cooler.27 

The  rate  coefficient  was  obtained  as  a  function  of  the 
electron  temperature  Te  by  integrating  the  energy  dependent 
cross  section  cr(Ecm)  over  the  electron  energy  distribution: 

8-777W,,  f°°  „  „  ' 

\  EcmO-(Ecm)e  Ecm'kTedEcm.  (9) 
(27rmekTe)  Jo 


B.  Branching  ratios  measurements 

The  following  DR  channels  of  Na+(D20)  are  energeti¬ 
cally  accessible  at  an  interaction  energy  of  0  eV: 


Na+(D20)  +  e~-+ 


NaOD+D+2.45  eV  (a) 
NaO+D2+ 1.79  eV  (/3) 
NaD+OD+O.96  eV  (y)’ 
Na+D2O+4.05  eV  (<5) 


(10) 


where  the  energies  correspond  to  particles  that  are  in  the 
ground  states  and  correspond  to  the  H  isotopes.  Isotopic 
shifts  are  small  compared  to  the  overall  exothermicities  and 
are  not  important  to  the  present  discussion,  so  have  not  been 
calculated.  It  is  interesting  to  note  that  the  recombination 
energy  of  this  ion  is  low  enough  so  that  there  are  no  exother¬ 
mic  channels  available  to  produce  three  neutrals. 

Since  neutral  particles  arising  from  the  same  recombina¬ 
tion  event  reached  the  detector  essentially  simultaneously, 
they  gave  rise  to  a  signal  proportional  to  the  full  beam  en¬ 
ergy  Eb,  which  made  it  impossible  to  distinguish  between 
events  originating  from' the  different  dissociation  channels. 
In  order  to  overcome  this  problem,  a  50  /mi  thick  metal  foil 
grid  with  70  fim  holes  was  inserted  in  front  of  the  detector. 
The  transmission  of  the  grid  is  P=  0.297±  0.015.  Since  some 
of  the  particles  were  stopped  by  this  grid,  the  spectrum  from 
the  MCA  showed  a  series  of  peaks  corresponding  to  the  re¬ 
combination  products.  For  instance,  the  products  from  each 
channel  could  be  either  detected  with  the  probability  P2  at 
the  full  beam  energy  or  with  the  probability  P(  1  -P)  at  the 
energy  carried  by  an  individual  particle.  Figure  4(a)  shows 
the  MCA  spectrum  recorded  at  a  collision  energy  of  nomi¬ 
nally  0  eV.  Fragments  originating  from  background  colli¬ 
sions  also  contributed  to  the  peaks  and  a  spectrum  was  also 
recorded  at  an  interaction  energy  of  1  eV,  where  recombina¬ 
tion  is  negligible,  see  Fig.  4(b).  At  this  energy  only  residual 
gas  collisions  are  considered  to  give  rise  to  a  signal.  There  is 
no  full  ion  beam  energy  signal  in  this  spectrum,  ensuring  that 
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FIG.  4.  (a)  The  signal  from 

Na+(D20)  +  e-,  measured  at  0  eV  in¬ 
teraction  energy,  (b)  Background  spec¬ 
trum  recorded  at  1  eV  interaction 
energy,  (c)  Background  subtracted 
spectrum  showing  only  DR  at  0  eV  in¬ 
teraction  energy.  Dashed  lines  repre¬ 
sent  Gaussian  fits,  corresponding  to 
the  different  fragments. 


electron  capture  from  the  residual  gas  is  negligible.  The 
background  subtracted  DR  spectrum  is  shown  in  Fig.  4(c) 
after  normalization  to  the  total  beam  intensity  during  each 
collection  run.  It  is  seen  that  the  Na  peak  is  broader  than  the 
D20  peak  [Fig.  4(c)],  even  though  the  corresponding  frag¬ 
ment  energies  are  similar.  However,  the  resolution  of  the 
D20  signal  is  related  to  the  energy  spreads  of  the  O  and  D 
signals  via  a  quadratic  sum,  and  from  the  fitting  [Figs.  4(b) 
and  4(c)]  it  can  be  concluded  that  the  width  of  the  D20  peak 
is  mostly  determined  by  the  O  contribution.  The  integrated 
intensities  of  the  peaks  are  similar. 

The  contribution  from  each  DR  channel,  Na,  Np,  Ny, 
and  Ns,  to  the  amount  of  the  observed  fragments,  ^d2o> 
^Na>  ^NaO  >  A NaOD .  and  /4NaE>2o ,  can  be  related  through  a  set 
of  linear  equations: 


^d2o  ' 

0 

0 

0 

p(\-py 

^Na 

0 

0 

0 

p(\-p) 

^NaO 

= 

0 

P(l-P) 

0 

0 

^  NaOD 

P(l-P) 

0 

0 

0 

_^NaD20_ 

A, 

p2 

p 2 

p2 

Na 


The  branching  fractions  na,  rip,  ny,  and  ns  were  deter¬ 
mined  after  normalization: 


In  order  to  determine  the  signal  levels  for  use  in  the  left  side 
of  Eq.  (11),  the  number  of  counts  of  each  fragment  was 
determined  by  fitting  the  data  to  Gaussian  curves  [Fig.  4(c)]. 

IV.  RESULTS  AND  DISCUSSION 
A.  Cross  section  and  rate  coefficient 

The  measured  cross  sections  for  DR  of  Na+(D20)  as  a 
function  of  center-of-mass  energy  ranging  from  about  1  meV 
to  0.1  eV  are  shown  in  Fig.  5,  where  error  bars  represent 
statistical  uncertainty.  The  systematic  uncertainty  corre¬ 
sponding  to  the  most  significant  errors,  viz.,  error  in  the  elec- 
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FIG.  5.  Cross  section  for  dissociative  recombination  of  Na+(D20)  as  a 
(12)  function  of  center-of-mass  energy  (solid  squares  with  error  bars).  Solid  line 
shows  the  best-fit  data  [Eq.  (13)]. 
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tron  cooler  length,  ion  and  electron  currents,  was  estimated 
to  be  20%.  A  least-square  power  law  fit  is  shown  as  a  solid 
line  in  Fig.  5.  The  equation  is 

o-(£cm)  =  (8.15±  1.63)X  1O~17EJbJ'37±0'01  cm2.  (13) 

The  cross  section  dependence  is  steeper  than  E~^  character¬ 
istic  of  the  direct  dissociative  recombination  mechanism.28 
The  same  tendency  has  been  observed  previously  for  other 
polyatomic  ions.29,30 

Rate  coefficients  for  thermal  distributions  of  electrons 
have  been  calculated  for  electron  temperatures  from  50  K  to 
1000  K,  which  can  be  fitted  by  the  equation 

I  j  \-0.95±0.01 

a=(2.3±0.32)X10"7  cmV1.  (14) 

The  observed  temperature  dependence  differs  from  the  T~0  5 
dependence  predicted  by  Bardsley31  when  only  the  direct  DR 
mechanism  is  taken  into  account.  This  indicates  that  DR  of 
the  Na+(D20)  ion  is  occurring  in  part  by  the  indirect  mecha¬ 
nism,  which  proceeds  through  vibrationally  excited  Rydbergs 
states  after  the  electron  has  been  captured  by  the  molecular 
ion. 

It  is  interesting  to  compare  the  present  results  to  those 
found  for  hydronium  ions,  H30+(H20)„.  Throughout  the 
DR  studies  of  water  clusters  with  an  H30+  core, 32-34  the  rate 
coefficients  have  been  measured  to  be  appreciably  higher 
than  found  in  the  present  study.  Estimates  of  the 
H30+(H20)„  DR  rate  constants  at  300  K  are  (0.5  +  2«) 
X  10"  6  cm3  s“ 1  for  n  from  1  to  4;  much  larger  than  those 
found  here.  For  the  bare  ion,  D30+,  which  has  the  same 
number  of  atoms  as  Na+(D20),  cross  sections  are  larger  by 
a  factor  of  2.5  to  3  than  found  here.  This  is  probably  related 
to  the  fact  that  recombination  energy  of  H30+  is  much  larger 
than  for  Na+(H20).  However,  the  electron  energy  depen¬ 
dence  is  almost  identical  at  E~lA6. 

B.  Branching  ratios 

The  following  branching  ratios  for  the  DR  of  Na+(D20) 
at  0  eV  interaction  energy  were  determined  by  solving  Eq. 
(11)  and  normalizing  the  results  Eq.  (12): 

«a(NaOD+  D)  =  0.05  ±0.01, 

«(8(NaO+D2)  =  0.03±0.01, 

«y(NaD+OD)  =  0,  (15) 

n^(Na+D20)  =  0.92±0.03. 

The  uncertainties  in  the  branching  ratios  indicate  the  varia¬ 
tion  of  the  result  when  the  uncertainties  in  the  grid  transmis¬ 
sion  and  the  fitting  of  the  peaks  were  included. 

The  DR  of  the  Na+(D20)  ion  is  a  two-body  dissociation 
process  tending  to  break  the  weak  Na+-D20  bond  (=92% 
of  the  total  reactivity)  which  as  shown  in  Fig.  2,  is  more  than 
twice  as  long  as  the  O-H  bonds.  A  DR  study  of  the 
D+(D20)2  cluster  found  similar  results,  where  the  cluster 
bonds  break  almost  exclusively.35  Channels  a  and  (3  produc¬ 
ing  a  D  atom  and  a  D2  molecule,  respectively,  require  one 
and  two  D-0  bonds  to  be  broken.  These  channels  account 
for  only  8%  of  all  DR  events,  indicating  that  the  covalent 


bonds  are  harder  to  break  than  the  Na+-D20  ionic  cluster 
bond.  Of  the  two  channels  where  covalent  bonds  are  broken, 
simple  bond  fission  is  70%  more  prevalent  than  the  channel 
that  breaks  two  bonds  and  forms  another.  The  D2  forming 
channel  may  involve  a  concerted  mechanism  where  both  of 
the  O-D  bonds  break  simultaneously  as  the  D2  bond  is 
formed  and  would  be  expected  to  be  enhanced  by  bending 
excitation  in  the  D20  molecule. 

The  channel  y  does  not  contribute  to  the  total  reactivity. 
The  confirmation  that  channel  y  is  absent  follows  from  the 
background  subtracted  pulse-height  spectrum  [Fig.  4(c)], 
where  contributions  to  NaD  and  OD  signals  are  negligible. 
Thus,  the  order  of  importance  of  the  dissociation  channels 
follows  simple  chemical  intuition.  The  descriptions  of  the 
chemical  nature  of  the  reaction  can  be  described  in  decreas¬ 
ing  order  of  importance  as  (1)  breaking  a  weak  ionic  cluster 
bond,  (2)  breaking  one  covalent  bond,  (3)  breaking  two 
bonds  and  making  another,  and  (4)  rearrangement. 


C.  Kinetics  of  the  chaperone  mechanism 

The  rate  expressions  for  Na+  and  Na+(H20)  in  the 
chaperone  mechanism  are 

4Na+] 

~~di  "  =  ~ Av[Na+3[H20][M]  +  *r[Na+(H20)][M], 

(16) 

4Na+(H20)] 

■  df  ■  =*/Na+][H20][Af|-Ar[Na+(H20)] 

X[M]-a[Na+(H20)][e"],  (17) 

where  kf  and  kr  are  the  forward  and  reverse  rate  coefficients 
for  reaction  (2),  respectively,  and  square  brackets  are  used  to 
denote  concentrations.  Two  simple  limits  can  be  described 
for  loss  of  Na+ .  The  first  is  if  the  effective  dissociation  rate 
constant  a[e~]  is  much  larger  than  the  effective  reverse  rate 
constant  kr[M\  and  the  second  is  vice  versa.  In  the  first  case 
the  rate  for  Na+  is  given  by  the  rate  of  clustering,  i.e.,  as 
soon  as  the  cluster  forms  it  recombines  with  electrons.  In  this 
limit, 

c/[Na+] 

r-j—  =  -kJLM}[  H20][Na+],  (18) 


where  kf  has  been  measured  in  He.15,36  For  the  atmospheric 
gases,  the  clustering  rate  is  expected  to  be  about  twice  that 
for  He  and  is  given  as37 


£/=  9.5X10" 301 


300 

\~t7 


1.5 


cm6  s"1. 


(19) 


For  the  second  limit  the  rate  expression  is  written  as 
4Na+] 


dt 


i^a[H20][<r][Na+]. 


(20) 


Keq  is  the  equilibrium  constant  for  Eq.  (2).  The  value  versus 
temperature  is 


Keq=136X\0~22l 


(*£_ 

A  H\ 
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cm 


(21) 
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where  R  is  the  molar  gas  constant,  and  AS  and  AH  are  the 
entropy  and  enthalpy  changes,3  which  are  deemed  to  be 
weak  functions  of  temperature.3,38 

With  the  above  expressions,  it  is  straightforward  to  cal¬ 
culate  lifetimes  for  atmospheric  conditions  with  respect  to 
each  of  the  three  processes  in  the  chaperone  mechanism. 
Since  the  atmosphere  is  cold  and  the  Na+-H20  bond  rela¬ 
tively  strong,  decomposition  can  be  ignored.  Using  standard 
atmospheric  conditions,  one  finds  clustering  times  on  the  or¬ 
der  of  (1-5)  X 105  s  at  80-85  km,  near  the  peak  of  the  Na+ 
layer.39  Higher  altitudes  have  even  longer  lifetimes.  In  this 
region,  the  electron  concentration  varies  between  ^  103  and 
105  e “/cm3.  This  translates  into  lifetimes  between  3000  and 
30  s,  respectively.  Comparing  the  two  lifetimes  shows  that 
clustering  is  rate  determining.  Since  pressure  and  water  con¬ 
centrations  increase  with  decreasing  altitude,  this  trend 
would  reverse  for  the  same  electron  concentrations  below  60 
km.  However,  in  the  natural  atmosphere  electrons  quickly 
form  negative  ions  at  that  altitude,  so  that  recombination 
would  be  mostly  ion-ion.40 

It  is  of  course  worth  comparing  the  chaperone  mecha¬ 
nism  to  three-body  recombination  involving  Na+.  There  are 
no  good  measurements  for  Na+  three-body  recombination. 
Therefore  we  estimate  the  rate  constant  to  be  that  found  in 
CaO  and  Johnsen,41  which  agreed  well  with  theoretical  cal¬ 
culations  from  Bates  and  Khare  42  At  200  K,  the  rate  constant 
is  ~2X  10“ 27 cm6  s“*.  Using  the  parameters  described 
above,  the  lifetime  at  85  km  is  3  X  107  s  for  an  electron  den¬ 
sity  of  105  cm-3.  Longer  lifetimes  apply  for  smaller  electron 
densities.  Thus,  while  the  rate  constant  is  larger  for  three- 
body  recombination,  the  greater  water  vapor  density  leads  to 
a  faster  loss  process  by  clustering  than  by  three-body  recom¬ 
bination. 

V.  CONCLUSIONS 

We  have  shown  that  the  DR  of  the  Na+(D20)  ion  is  an 
efficient  process.  The  cross  section  was  measured  over  col¬ 
lision  energies  ranging  from  1  meV  to  0.5  eV  and  shows  a 
E~137  dependence.  The  DR  rate  coefficient  was  deduced  to 
follow  (2.3±0.32)X  10_7(7’g/300)”°'95±0'01  cm3  s_i,  indi¬ 
cating  the  presence  of  the  indirect  mechanism.  Rate  coeffi¬ 
cients  for  a  chaperone  mechanism  have  been  calculated  for 
limiting  values.  The  complete  branching  ratios  were  deter¬ 
mined  at  the  minimum  center-of-mass  energy.  The 
Na+(D20)  ions  dissociate  only  into  two  neutral  fragments 
and  Na+D20  channel  is  the  most  dominant  (92%). 
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